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ABSTRACT 

Exact image method, recently introduced for the solution of electromagnetic field prob- 
lems involving sources above a planar interface of two homogeneous media, is shown to be 
valid also for sources located in complex space, which makes its application possible for 
Gaussian beam analysis. It is demonstrated that the Goos-Hanchen shift and the angular 
shift of a TE polarized beam are correctly given as asymptotic results by the exact reflec- 
tion image theory. Also, the apparent image location giving the correct Gaussian beam 
transmitted through the interface is obtained as another asymptotic check. The present 
theory makes it possible to calculate the exact coupling from the Gaussian beam to the 
reflected and refracted beams as well as to the surface wave. 
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INTRODUCTION 


The Gaussian beam is an important model for radiation from large aperture antennas 
and laser sources, from which the energy is sent in a relatively narrow space angle. Intro- 
duction of a complex space dipole to mathematically represent a source for the Gaussian 
beam by Deschamps [l] made it possible to extend analytical results, derived earlier for 
real space sources, to Gaussian beam excitations. The basic Gaussian beam originates 
from a dipole in complex space, while multipoles of higher order can be shown to produce 
more complex Hermite-Gaussian and Laguerre-Gaussian beams [2], [3], [4). 

The basic problem in Gaussian beam analysis is the problem of reflection and refraction 
of the beam at a planar interface of two homogeneous dielectric media, where losses might 
be present. This problem has been treated both through complex space dipole method 
[5] and other methods [6], (7], [8], The general field expressions in the exact formulations 
are involved, but asymptotic considerations for narrow beams allow simpler analysis of 
reflection and transmission problems and lead to the well-known shift phenomena of the 
reflecting beam. The most famous of these was first demonstrated by Goos and Hanchen 
[9] in 1947 and analyzed by Artman [10] in 1948, involving a parallel shift of the beam 
reflected from an interface in denser of two lossless dielectric media. Other shifts which can 
be defined for narrow beams are the angular shift [5], [11] and the focal shift [12]. Infinities 
arising in the idealized theory for the Goos-Hanchen shift can be avoided through a more 
realistic analysis [6], [13], [17], 

The purpose of the present paper is to give a more general account of the exact image 
theory, introduced elsewhere for sources in the air above a planar interface of a homoge- 
neous medium [14], [15] (the Sommerfeld problem) and allow sources to be in complex 
space. Thus, the theory is applicable to global Gaussian beam analysis in exact form. In 
comparison with other exact methods applying Sommerfeld integrals, the present theory, 
dealing with sources and their integration, gives more physical insight and does not resort 
to special integration procedures dependent on the field point. For the application of the 
exact image theory, functions characterizing the image sources are needed and methods 
for their computation have been presented in [14], [15], These functions need only be cal- 
culated once in the computer memory, after which the field calculation involvs converging 
well-b haved integrals. The exact image of a point source is a line source, which is nor- 
mally in complex space. As checks to the theory, known asymptotic expressions for the 
reflected beam shifts and the apparent transmission image location are shown to arise as 
special cases when the exact line image is approximated by a point image. 


REFLECTION IMAGE THEORY 

The notation applied here is based on Ref. [14] except that we consider two half spaces 
with respective parameters miMO? ejeo for u • t > 0 and M2M0> £ 2 e 0 f° r « • r < 0. In [14], 
[15], the medium u -f > 0 was assumed to be air, but generalization for any medium pair is 
straightforward, if ins' d of k we write fcj — ky/fi jej, where k — w^/Jiqc o, and everywhere 
understand that /i — t = C2/ e l- To avoid unnecessary, although not excessive, 

complication, the theory here is limited to dielectric media only with Ml = M2 — !• 

As another generalization, the original source is taken to be in complex space. Because 
the theory in [14], [15] was analytically derived and no use of the tacit assumption of real 
source location was made, the final expressions can also be applied for complex source 
locations. The same idea has been applied earlier with success when field expressions, de- 
rived for problems with sources in real space, have been generalized for sources in complex 
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space (1 6], resulting in solutions for problems with Gaussian beam excitation. 

Let us consider, for simplicity, the source J(r ) ~ vlL6(f - (uh +jb)), which is a dipole 
with the direction of the unit vector v and location at the real height h from the interface 
u • r ~ 0 and imaginary depth jb. When calculating the fields in the half space 1, the half 
space 2 can be replaced by the exact image source, as was shown in {14]: 


./,(?•, p) = /i(ppri(r)- (/f(/>) + (1) 

Here p is an integration variable and the image function / f (p) can most conveniently be 
calculated through the following Bessel function series m 



Further, 6+(p ) denotes 6(p — 0+), t a component transverse to u, and c the reflection 
operation 

a c = C’ • o, Jc(f) = X' • J(C • r), (7 = 7 — 2uu. (3) 

The field from the image source can be written as a fourfold integral over space and the 
parameter p: 

E(f) = -ju>p.Q j f G\{D) • Ji{f' ,p)dV'dp, (4) 

JV J 0 

with 

= /= 1 \ e~i k ' D 

D(r,r',p) = y/(f - f' + up/jB) ■ (r ~ r' + up/jB), B - ky/e 2 - ej. 

The image current expression contains Bessel functions J n (p ), which are convergent only if 
the integration parameter p is real. For a point source, (1) can also be written so that the 
integration parameter p is absent, because in the field integral (4), the coordinate : — u-f 
and p are related through the distance function D in the argument of the Green function. 
The expression ( 1 ) can thus be written 

Jj{f) - ._/£-——.{}(£ . r)6 + (r + uh - jb c ) - jk\y/(. - II L [t>c/i(p(z))$(P - ' jb)~ 

-u(u • v)ff(p(s))6(p - jbt) - ttp-i-i/f(p(i))tt • ^(uc * V)(5(r 4- uh - >6 C )| , (6) 

with 

p(z) = -jB(z + h+ju ■ b). it (7) 

To obtain an exponentially converging Green function, we must select the branch of the 
distance function D so that Irn[k\D ] < 0. Also, to obtain a converging //image current 


function, the path of z integration must be chosen so that p in (7) is real and positive, in 
which case all the Bessel functions in (2) converge. This means that the image line must 
start at the point r = -uh + b c in complex space and lie parallel to the complex z plane 
with 

arg(z -f h -ft ju-b) = arg(jB). (8) 

Requiring that the additional condition Rcjs] ^ 0 be valid, the complex distance function 
satisfies D ^ 0 for all field points in Re\z) > 0 and the field integrand in (4) is nonsinguiar 
[14]. This condition defines the branch of B } or the square root y/e 2 — cj through 

Re[jy/e 2 - cl] == Rc[v/ti' — t'j] > 0, (9) 

i mplyin g for example, that if cj ,<2 ~ f«i <1 «i are real, then v/c 1 = ~js/\ - e, provided 
y/l - e > 0. For some combination of lossy medium parameters there may be doubt of 
choosing the correct branch of the distance function D to obtain the best convergence [18]. 
A close study of the integration path on the complex ~ plane will reveal that a branch cut 
line, starting from the branch points at z ' ■ == s ± jp, defining the converging branch of the 
Green function, may be crossed for some combination of parameter values ej, ( 2 , when p 
moves from 0 to 00 . In this case it is possible to take the image current line in the ’wrong’ 
half space with Re\z') > 0, which makes the integrand converging again. This question is 
the subject of a forthcoming paper. 


REFLECTED GAUSSIAN BEAM 

There does not seem to exist a global definition for the Gaussian beam field, instead, 
the term ’’Gaussian beam” is understood as an asymptotic property of radiation fields 
close to the axis of the main radiation. Thus, many fields with the same asymptotic 
quadratic exponential behavior are called Gaussian beams. One example is the field from 
a point source in complex space with the imaginary part of the position vector, 6 = 
ubcqsd — UybsinS. 

The radiation beam of the point source is obtained in the direction of —b. The reflected 
beam is obtained from the image source (6), and, for points far enough from the interface, 
the field can be approximated by a Gaussian beam, if the image source can be approximated 
by a point source. To demostrate the validity of the exact image theory in this generalized 
sense, we show that the Goos-Hanchen shift as well as the angular shift of the reflected 
beam are obtained as asymptotic results from the exact image field expression. To keep 
notation simple, let us consider a two-dimensional problem with the line source J(f) — 
u r I6(y + jbsin0)6(z — h — jbcosd), where 0 is the angle between u and b. This assumption 
simplifies the image source (1) into 

ji{r,p) - fl(p)Mr) = fi{p)vxI6(y + jbsin0)6(z + h + jbcos0). (10) 

The Green function (5) integrated in x direction produces the two-dimensional Green 
function in the medium 1, 

5i(D) « (7+ pVV)0,(D), G,(fl) = ijr< 2) (*,X>), 


( 11 ) 
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( 2 k 

where //q (j) is the Hankel function. 

The field calculated from the exact image source is valid in any point in the upper half 
space, but since we are interested in Gaussian beam properties, let us consider the far field 
with |At|Dj ! >> 1 , in which case the image can be approximated by a point source. Also 
the asymptotic expression for the Hankel function can now be used: 


HfkkiD) « 



-ihD 


( 12 ) 


D here represents the com plex distance between the field point f and the complex inte- 
gration point r* + jup/h'i y/e - 1 , where r' = ~uh + jb c , if p integration is done separately. 
Because the image function /j(p) is decaying, the effective p integration range can be 
regarded as small with respect to |r — f'| and we can write 


n n JPQ n fT- -n -«V « 1 _ V - / 

Z> » A> - T -'i Do = y/{r - r 1 ) • (r - r 1 ), q = - 7 ==, w-~=—. 

kl y/( - 1 Do 


(13) 


Applying (12), the field integral can be written on the plane * = 0 as 


/l(p)c“ M dp. (14) 

If (14) can be written in the form 

E(r) -ju>p 0 v x IG j ( Do * ~ s c " t * oi [°° /i(p)e" w °dp, (15) 

J 0 

in the vicinity of the radiating direction w xz w 0 , where w 0 corresponds to a field point f 0 
on the axis of the reflected beam, the field can be thought of as arising from the current line 
shifted by the vector 5 from the mirror image location r\ Approximation (15) is obviously 
possible, if (14) can be expanded as a power series in terms of the small difference vector 
w — w 0 - This is possible only for narrow enough beamwidths. After some steps of Taylor 
expansion, the expression for s can be written in the form 


E(f) RS -juipo^xIG 


h (Do) f 
Jo 


Jo* pfi(p)e~ M °dp 
1 /o°° /i(p)« _w °dp ’ 


(16) 


where qo — u * u>o/ y/e ~ 1. To obtain an expression for s, the following integral identities 
are needed: 



/l(p)e ™dp 


q - y/g 2 + 1 
q+ yjq 2 + 1 ’ 


(IT) 


/ Pfl(p)e W dp = 
Jo 


2 q-\/q 2 + 1 

y/q 2 + lq + y/ q 2 + 1 


(18) 


Thus, we have 


s — 


2 ju 


2 ju 


2 ju 


kiy/7=iy/q 2 + l W * 2 “ <1 + «l(fi • ®o) J *|V« — sin? 


\ 


(19) 


5 


with the branch of the square root so defined that "~ain^0\ < 0, {14], (19) explaiiis^ 

the Goos-Hanchen shift and the angular shift when the Ta ylor ap proximation condition is 
valid. In fact, because w a = ~b c /b , we have qo — coa0/ y/T~ T, and (17) represents the 
TE reflection coefficient of a plane wave coming at the angle 0. Let us consider the two 
possibilities with real c: 

1) c < ain“0, which can only happen for cj > fo. In this case, 5 is real and thus the 
real part of the location vector is shifted by 


a 


kiVsin^O — c’ 


( 20 ) 


which means a parallel shift of the reflecting beani, Fig. 1. This is the well-known Goos- 
Hanchen shift. Because 5 is in — y direction, the parallel shift is ssinS. 

2) e > ain“0. In this case, the square root is real and the shift 5 is imaginary = jau, 
which means that, the image line is shifted an imaginary distance in the u direction from 
the original location. Because the imaginary part of the position vector determines the 
direction of the beam, the beam is shifted angularly from its mirror image direction w 0 to 
the direction determined by b c - js, Fig. 2. If the shift angle A0 is small, it satisfies 


tar i(A0) =5 — • £cos0 = — - : — ~^=====. (21) 

b k x bsft - sin?0 

The previous expressions are only valid for sufficiently narrow beams, which presumes 
sufficiently large values for b, They are not valid if the Taylor expansion is not applic able, 

which happens at and close to the branch point of the square root function y/?o + 1, i.e., 

at sin“0 « sin 2 0 c — c, where (19) would predict an infinite shift. This is the definition of 

the critical angle 0 C , for which yjq^-r^ = 0. For angles 0 0 C W« can write from (14), 

< 17 ) ■ ’ 

E(f) w -ju>/i 0 fi*/G 1 (I)o)t" 2j V 2 <^-^>/^ r -(22) 

which cannot be described by a simple shift of image source, because the reflecting beam 
is distorted and not exactly of Gaussian form. To find out the shift of the maximum of 
the reflecting beam, some numerical analysis must be made, as in [6], [13], [17]. 

For complex c values, there are both real and imaginary parts for the vector 3, which 
means a combined Goos-Hanchen and angular shift. In this case, the beam does not, enter 
at the critical angle, because it is now complex. For beams narrow enough, the real and 
imaginary parts of 3 of (19) determine the Goos-Hanchen and angular shifts, respectively. 

This simple test demostrates the applicability of the exact image theory for the Gaus- 
sian beam analysis. The field of the reflected beam can be calculated at any point from the 
exact formulation (4). This also includes the coupling to the surface wave, which however 
is small when 0 is not ss 0. 


TRANSMISSION IMAGE THEORY 

The image theory for fields transmitted through a planar interface was developed in 
ref. [15] and the generalization discussed for reflection image theory also applies here. To 
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express the result, the upper dielectric half space is replaced by the transmission image 
source: 


srM. 3(f) (“<+« + + J,(f)(/i + (p) + fJ(fl 2 , P )) + 


«[F,( B-.,p) - F,{B:,p)]U‘U,p)T, • J( r). (n) 

Here, the original source is J(r) and z denotes its coordinate. Further, we have 

B — k\y/t~ l - ky/t2 ~ 7i, (24) 


H(z } p)~Jz* + (p/B)'>, 


F ( (t,p) s ^-Jo(p) 


__ V s / e - 1 1 - V(p/t ) 2 + l \ 

1 V + U + ^p/r) 2 + 1/ 


*^2m(p)t 


and 1 in functions H and F denotes differentiation with respect to p, 

The field is ohu^ted from an integral similar to (4) but the medium is now to. Also, 
the distance function D is here defined to be 


D(f,r',p) = ~ - uH(z',p)) • ( f-p ' - uH(z\p)). 


Again, to obtain a converging image function, the parameter p must be real. To obtain 
diverging Green function, Im^D] must be nonpositive, which defines the branch of 


a converging Li reen lunction, lm\k 2 U\ must be nonpositive, which dehnes the branch ol 
the D function. Also in this case, for certain lossy media, the branch cut of the Green 
function may be crossed if the image is restricted to the half space JSe^] > 0, leading to 
nonconverging Green function. This problem can be overcome by taking the other branch 
of the H(z,p) function. 


TRANSMITTED GAUSSIAN BEAM 


Without delving more into the theory itself, which can be found in [15], let us apply 
it to the same line source as in the previous Section. The transmission image source can 
be written in the form 


Ji(f, P ) = ux/(Mp) + *i(fl/t\p))% + jbsin0)6(z - /»'), 
where the function F{ is defined as 


Fi(Bh',p) 


+ ft') 


((H +h'fj,(p) + (jf 


H s H(h‘,p) = Jh* + (p/fl)2, ft f = ft + jbcosB. 
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The convergence condition Im[p\ ~ 0 defines a path : = JH(p), which consists of a part 
of hyperbolic curve in the complex z plane [15]. To obtain an asymptotic result from the 
exact image field solution, let us once again study the far field of the beam, Making use 
of (12), the following approximation for |f| large must now be applied instead of (13), for 

x » 0, s < 0: 

D *s Dq + Do - \f(y+jbtinB)- + q *= \z\/Do coaB 2 > (31) 

Thus, the field in medium 2 can be written in the approximative form 

E(f) * ~ju>p 0 u r JG 2 (Do) (l + /°° F' 1 (Bh' } p)e~ ik * qH Mdp\ . (32) 

Here, the Green function is defined similarly as in (5), for the mediuny 2 with k% 
replacing k j. Applying an integral identity [15, eq. (17)], we can write for the expression 
in braces in (32): 


1 + 


rniBk', P »-»*v<p 

Jo 




7 + v/ 7 2 ~ 

Use of this and (see Fig. 3 for the different parameters) 

Dq r + jbainBainB-i fs 4- (fainBo, 

with 




<f = d + jbainB, 

makes it possible to write (32) in the following form: 


(33) 


(34) 

(35) 


1 -r 2A'2 Cos^2 

£(r) = , 


-jhrt e ~jk 1 linfi > d' e -jk t coa0 1 h! 


Here, we have made the far field assumption r 2 > h,d,b, otherwise the expressions would 
have been much more complicated. (36) can be interpreted in terms of geometrical optics 
rays of the Gaussian beam, with a phase shift in each medium and a transmission coefficient 
at the interface. From this, an equivalent source point location can be defined from the 
divergence of two adjacent rays close to the beam axis B\ = B. Taking the combination 
of the last two exponentials in (36), we may require, that the corresponding expression 
for the image case has the same change when the angle B\ differs from the axial angle B. 
Writing for the exponential expressions 

= ~jki{d' ainBj + h'cosBj) — —j —~ — •+• Aj5cos(0i — B), (37) 

coao j 


. M 2 


^2 = — jA'2(d!>sin^2 + hoCoaB^) — —j — + A?2^2 cos (^2 ~ ^o)» 

cosv 2 

k\ sinB — k2i>inBQ, 


(38) 

(39) 


we may require that phase and amplitude changes along the plane r = 0 be the same for 
the original and the image case. This implies that the differentials of and $2 are the 
same: 


( 40 ) 

When supplemented bythe condition between the differentials d0 j, d$ 2 , arising from Snell’s 
law = A'2 ^i»^ 2) (40) gives us the following expression for the unknown ;h % : 


h h^l } al e, 


0 . 


k\coa*0\ 

Because of L\ *= hjcoa0\ , £>2 ™ hi/cotBi, the corresponding result for £2 * s 


(41) 



j ^2 coj2 ^2 
1 A^COS 2 #! ’ 


(42) 


which was also given in (5). To obtain an expression for 62, defining the apparent source 
location in complex space, we must study the real parts of (37), (38). The direction of the 
vector is known from the main beam direction in the half space 2. The magnitude 62 is 
not obtained from (40), since the first differentials vanish identically at 0\ = 0. Taking the 
second differentials of and ¥2 and equating the real parts gives us the result 


62 
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A’2C£»i^^2 

AjjCOS 2 #! ’ 


for 0\ ss 0 , 


(43) 


which condition was also given in [5]. 

Thus, the location of the approximate image source can be obtained through two differ- 
ent approaches: the present exact image theory and the saddle-point asymptotic analysis 
of Fourier integal representation of the transmitted field in |5]. The present approach has 
the advantage of working with sources, which gives a more physical insight for the problem 
in addition to being exact. 


CONCLUSION 

In this paper, the exact image method, previously introduced for problems involving 
sources in real space with two homogeneous media and a planar interface, is extended 
to problems involving sources in complex space. This is possible because the analytic 
functions of real space source position can be extended to functions of complex space 
source position, which idea has been applied before with success in diffraction problems. 
The complex source point theory gives a possibility to analyze Gaussian beam problems 
with a simple source instead of a complicated source in real space. The method has been 
tested in this paper by analyzing asymptotic (far held) expressions for a reflected and a 
transmitted Gaussian beam, resulting in well known Goos-Hanchen and angular shifts for 
the reflected beam and apparent position for the transmitted beam. The present method 
suggests itself for more exact calculation of the fields in these beams. 
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Figure K, 

Goos-Hanchen shift of the Gaussian beam arising from an approximate image point 
source with a real shift vector S from the mirror image point *~uh 4- jb c . 



Figure li. 

Angular shift of the Gaussian beam arising from an approximate image point source 
with fin imaginary shift vector s from the mirror image pointVfih + J*o 
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